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Abstract 

Properties of a given symmetry group G are very important in investigation 
of a physical system invariant under its action. In the case of finite spin systems 
(magnetic rings, some planar macromolecules) the symmetry group is isomorphic 
with the dihedral group Djy. In this paper group-theoretical 'parameters' of such 
groups are determined, especially decompositions of transitive representations into 
irreducible ones and double cosets. These results are necessary to construct matrix 
elements of any operator commuting with G in an efficient way. The approach 
proposed can be usefull in many branches of physics, but here it is applied to finite 
spin systems, which serve as models for mesoscopic magnets. 
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1 Introduction 



One of standard ways of quantum magnets modeling is to consider a finite set of nodes 
carrying spins Sj, 1 < j < N, with a spin number s. At present such models are also used 
to investigate properties of magnetic macromolecules as Mni2acetate |l|]. The so-called 
Ising states or configurations, i.e. tensor products 



of one-spin states \mj) form a natural basis B of the quantum space of states L with 
the dimension dimL = (2s + 1)^. In a general case, however, these states are not 
eigenvectors of the Heisenberg-like Hamiltonians. Therefore, it is necessary to solve an 
eigenproblem with dimL equations what, even using contemporary computers, is almost 
impossible for ^ 30 and s > 1. This is especially important when one wants to 
investigate spin correlations in the ground state and some of the first excited states of 
antiferromagnetic systems, since many methods allow to determine eigenvalues, but it 
is much more difficult to find eigenvectors. Among others, it is caused by the numerical 
instability of necessary algorithms in the case of a very large number of linear equations . 
The problem is partially solved by taking into account operators commuting with a given 
Hamiltonian: e.g. the 2;-th component of the total spin and irreducible representations 
of the Hamiltonian symmetry group 0, If the spin interactions are isotropic, then also 
the total spin number may be taken into account. However one has to find solutions of a 
set of homogeneous linear equations arising from the eigenproblem for (in finite spin 
system eigenvalues of are known). In many interesting cases, for example investigating 
the ground state of bipartite antiferromagnets, these solutions can be found exactly if one 
is able to cope with very large integers. Such tools are provided by multiple (or arbitrary) 
precision libraries like the GNU Multiple Precision Arithmetic Library (GMP or GNU 
MP in short) written and developed by T. Granlund However, such procedures work 
efficiently if the matrix of an operator (the squared total spin or the Hamiltonian Ti 
in the presented considerations) are written in the symmetry adapted basis. 

The considerations presented are illustrated by examples related with finite spin mod- 
els, but a similar method can be applied in any case when: 

• the space of states L is spanned over orthonormal vectors (0), i.e. it is the A^th 
tensor power of the n-dimensional one-particle space; 

• the considered operator Ti, acting in L, commutes with a given permutation repre- 
sentation P of a group G (the Hamiltonian symmetry group). 

Let a given group G permutes vectors of a (finite) basis B of the quantum state space 
L. In the previous paper general formulas for vectors of the irreducible basis and 
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Figure 1: Numbering of ring nodes and elements of the group Dg (the six-fold axis Cq is 
denoted by the full hexagon) 

matrix elements of any operator commuting with all permutation operators P{g), g E G, 
have been presented. The aim of the present work is to investigate the dihedral groups 
Djv and obtain general analytical expressions in this case. These groups have been chosen, 
since they describe symmetry of molecular rings with identical spins s In fact, these 
rings are also invariant under the reflection ah in the horizontal plane, so the groups D^vh 
should be considered. These groups are important, for example, constructing symmetry 
coordinates of vibrating molecules or crystals, but in the case of spin models the action 
of ah is trivial. Considering planar molecules one should introduce groups Cnv, which 
are elementary 'bricks' for finite space groups 0. However, groups Dn and Cnv are 
isomorphic and the abstract group D^v is considered as a representative of isomorphic 
groups with the following generation relations 

D^=(C^,?7o), = Ul = {&M^ = E, 0<j<iV. (2) 

The dihedral group Dn consists of 2N elements: powers of the A^-fold rotation Cn 
and perpendicular to it the twofold rotations Uj = Cj^Uo, < j < N. It is easy to write 
down general formulas for multiplication (all calculations are done mod A^) : 

^N^N ~ 1 Cj^Uk = Ujj^k , UjCj^ = Uj-k , UjUk = C'l^ . (3) 

The labeling scheme for = 6 is presented in Fig. |I[ 

The dihedral groups are considered in Sec. H; there are many differences between cases 
of A^ odd and even, so the results are presented separately. In many actual applications 
one deals with bipartite antiferromagnets so the case of A^ even seems to be more 
important, however the other case (A^ odd) is discussed for three, at least, reasons: (i) the 
ground state energy of infinite linear antiferromagnets is approximated from below and 



above by those obtained in the cases of even- and odd-sized rings, respectively (ii) 
there are still some open problem related with linear ferromagnets {e.g. critical exponents 
near T = OK |^); in this case the parity of is (almost) irrelevant; (iii) the presented 
considerations and results are more complete. In Sec. ^ a partial solution for six spins 
s = 1 (a ring invariant under Dg) is discussed. The paper ends with some conclusions. 

2 Dihedral groups Dat 

Let us recall the formulas obtained in the previous paper |^ to check what group properties 
should be determined. The irreducible basis consists of vectors \U fi; Tr'y), where: 

• f/ C G is a subgroup of G representing a conjugacy class f/; 

• H is a representative of an orbit of type U ; 

• r is an irreducible representation (an irrep) of G, the repetition index r distinguishes 
copies of r and 7 = 1, 2, . . . , [F] labels vectors of T. 

The orbit G{fi) can be identified more precisely if one introduces a label denoting a non- 
ordered partition [k] such that fi G B[k]. Note that [k] identifies also the magnetization 
M. However, the formula for vectors of the irreducible basis does not depend on [k] and 
is given as 

\G/U\ 

\u^^■,Tn)= y: &r^(f/)|t//x^?,), (4) 

where the elements gj, 1 < j < \G/U\ with gi = E, represents the left cosets gU G G/U. 
The coefficients , determining a unitary transformation, can be expressed by matrix 
elements of the irrep F and the reduction coefficients Al^^: 

hY\u) = i^^T E r;,(^?.)«0* • (5) 



/3=1 



Matrix elements of any operator, e.g. a Hamiltonian ?i, commuting with all P{g) are 
grouped in blocks labeled by F and depend only on [/, /z and r: 

/|rr|\l/2 \G/U\ 

hu,ry.s{T) = ^ {UMmVi^E)B^r{9j) , (6) 



\y\j 



where 

i?.^/"(^.)=E(^rj)Tc.,(<7,x'- (7) 

It is helpful to consider non-ordered partitions [k], because matrix elements are equal to 
zero if partitions [k] and [k'] are not joined by an edge in the graph introduced in Sec. 5 



of the previous paper 0. Introducing the notion of double cosets VgU this formula can 
be modified and written as: 

\V\G/U\ 

huf^ryusir) = {\U\\V\)-'/' \Vg,U\Bg'''^{g,){U^igd\n\VuE), (8) 

d=l 

where elements gd represent double cosets VgaU, or orbits of the natural action v{G/U). 

All structures appearing in these formulas (cosets, irreps, reduction coefficients etc.) 
in the case G = D^v are discussed in the following subsections. 

2.1 Classes and irreps 

The multiplication rules (^ in the dihedral group Dn lead immediately to classes of 
conjugated elements: 

N odd : {E}, {C^, J = 1, 2, ... , , {f/o, t/i, . . . , Un-i} ; 

AT even: {E}, {C^, C^'}, j = 1, 2, . . . , ^ - 1 , 

{C2 = C^'^], {f/o, . . . , Un-2] {f/i, f/s, . . . , Un-i). 

Therefore, there are (A^ + 3)/2 classes (and irreps) for A^ odd and A^/2 + 3 for A^ even. In 
the special cases, A^ = 1,2, the groups Dj^ are Abelian. 

In both cases, A^ odd and even, there are two one-dimensional irreps, denoted Ai (the 
unit irrep) and A2, with characters 

X^H^?) = 1, V(7GD^; (9) 
X^HC^) = 1, X^\U,) = -l, J = 0,1,..., AT -1. (10) 

For A^ even there are two more one-dimensional irreps, denoted Bi and -B2, with characters 
(j=0,l,...,A^-l) 

x^nc^) = (-iy, x^Hf/.) = (-ly ; (11) 

X^'KC'n) = {-ly . X^^(f/.) = -(-l)^ (12) 

Note that the element C2 = C^^^ behaves in different ways for A^ = AM and A^ = 4M + 2. 
In the first case one obtains X'^^(C*2) = X^^{C2) = 1, whereas in the second case this 
character equals —1. 

The other irreps, denoted Ei, I = 1,2,..., A^^;, are two-dimensional; there are A"^; = 
(A^— 1)/2 and A^^; = N/2 — 1 such irreps for A^ odd and even, respectively. Their characters 
are given by the same formulas 



X^' (C^) = 2 cos(27rj7/Ar) , x""' m = 



(13) 



To complete the presentation, matrices of the representations Ei in the standard basis 
{\Ex), \Ey)} for the generators Cjv and Uq are written down: 

2.2 Subgroups and conjugated subgroups 

Let be written as a product of r different primes pf. 

i=l 

Then there are I> = ni=i(Q;i + 1) different divisors k of N, each of them can be written 
as {k — 1 corresponds to all Pi — 0) 

r 

i=l 

The first family of subgroups consists of the cyclic groups C^, where /c is a divisor of 
A^. These groups are generated by elements C^, where k = N/k is the so-called co-divisor 
of k. For k = 1 one obtains the group Ci = {E}, which has no generators; in this case 
K — N. Each subgroup contains k elements in the form Cf^ = C^,p = 0,l,...,k — 1. 
Since these groups are normal divisors of Djv, then each is conjugated to itself only. 

The second family of subgroups is formed by the dihedral groups D^. However, for 
each k, i.e. for each generator C^, one has freedom in a choice of two-fold generator 
Uj. The subgroups isomorphic to are generated by and one of the rotations 
Uo,Ui,. . . , Uk-i. These subgroups will be denoted hereafter as 

^k = {CN,U,), q = 0,l,...,n-l. 

Each of these groups contains, except for the rotations C^^, twofold axes Uq+p^, p — 
0, 1, . . . , /c — 1. In the case of odd k, these elements belong to the same conjugacy class, so 
all subgroups are conjugated. This conjugacy class is represented by Dj. = {C^, Uq). 
When K is even, then all two-fold Uj rotations in the subgroup have indices with 
the same parity. Therefore, there are two classes of conjugated subgroups: the first is 
represented by = {C^, Uq) and the second by = {C%, Ui). If N is even then it can 
be written as 2°'N' , where N' is odd. In the same way all divisors k ol N can be written 
as 2^k'. II (3 = a then corresponding k = N'/k' is odd and, therefore, there is only one 
class of conjugated subgroups D^. Two separated classes, and D^, are obtained in the 
other cases, i.e. when (3 < a and k = 2°'~^N' /k'. Let V = T>/{a + 1) be the number of 
divisors of N' {V = V ii N is odd) then: 

• there are V classes of dihedral subgroups in type D^; 




Figure 2: Lattice of conjugated subgroups for the dihedral group Dg; hues correspond to 
the group-subgroup relations 

• in each type D° and there are oD' classes of dihedral subgroups. 

Therefore, there are {a+l)V' cyclic subgroups and {2a+l)V' conjugacy classes of dihedral 
subgroups. The lattice of conjugated subgroups for the group Dg is presented in Fig. |[ 

The theorem presented in Sec. 2.3 of the previous paper immediately excludes 
some classes of subgroups if vectors of the considered basis are mappings, as in the case 
of magnetic configurations. The first excluded subgroup is the cyclic group Cat, since 
there is only one orbit X = {0, 1, . . . , — 1}, so Yij = and P(J)n) fl Sx = 
P(J)n) 7^ P{Cn)- In the case of even two other subgroups are excluded: the dihedral 
subgroup Djy^2 and the cyclic subgroup CAr/2. In the first case only one orbit X is again 
obtained, whereas CAr/2 decomposes X into subsets Xq and Xi containing nodes with 
even and odd indices, respectively. Then one obtains P(I)n) fl (Sxq x = P(D%/2) 
P{Cn/2)- Therefore, these subgroups can be omitted in the further considerations of spin 
configurations. For example, when = 6 subgroups Cg, Dg, and C3 are excluded. 

2.3 Cosets 

For each representative f/ of a class of conjugated subgroups U left coset representatives 
have to be chosen and fixed. The index |DAr|/|D;| = A and the rotations C^r, j = 
0, 1, . . . , A — 1 will play this role. In the case of cyclic subgroups one needs additional 
A representatives. They can be chosen as Uj, j = 0, 1, . . . , A — 1. The choice of coset 
representatives is fixed in the further considerations. 

Considering double cosets V\Di^f/U we determine at first their number in different 
cases using results presented in Sec. 3 of the previous work 0. If both U and V are cyclic 



even 



then only classes containing rotations Cj^ have to be taken into account and one obtains 

|C,\D^/Q| =2gcd(/«,A) =2x, 

where 'gcd' denotes the greatest common divisor and n = N/k, A = N/l are co-divisors of 
k and /, respectively. Note that for mutually prime co-divisors k, A there are two double 
cosets only: one contains rotations C^, whereas the second — elements Uj. If \^ = 
(or U = Di and V = Ck), then the only change is the order of one subgroup, so the result 
is 

|CADjv/D,| = |DADjv/Q| = x. 

If both subgroups are dihedral then the cases of odd and even have to be considered 
separately. Moreover, in the case of even N a dihedral subgroup can be in one of three 
types: D^, or D^. In the first case it contains elements Uj with even indices j, in the 
second case all j's are odd and in the last case odd and even indices j are present. The 
final results follows: 

• iVodd 

|D,\D^/Dz| = (x + l)/2; 

\Bl\B^/D^\ = \Dl\D^/Dj\ = (x + 2)/2; 

\Bl\Dj,/Dj\ = x/2; 
\Dl\D^/Di\ = \Dl\DN/Bi\ = \Dk\BN/Bi\ = (x+l)/2. 

The above presented formulas are symmetric, i.e. the subgroups U and V can be 
transposed. However, in a general case a set of double cosets V\G/U is not symmetric, 
i.e. it differs from the set U\G/V (they have the same orders only). The same can be 
said about representatives Qd of double cosets. Therefore, both cases should be considered 
separately and for each pair V, U we obtain two sets of double coset representatives. On 
the other hand, the operators considered are Hermitian, so hu^r,Vus = ^yvs,(7/xr ^^^d, 
therefore, it is enough to determine one element of each pair. Nevertheless, for the sake 
of completeness all cases are considered below. 

Since the cyclic subgroups are self-conjugated, then for a given V all double cosets 
VgCi have the same order (see Eq. (13) in 0): 

IWI = |C,<,V|=|^. 

The stabilizer in the denominator for V = or Cfc is equal to Cgcd{k,i)- Therefore, 

kl 

\CkgBi\ = \DkgCi\ = 2\CkgCi\ = 2 = 21cm(A;,/) , 

gcai^fc, L) 



where '1cm' denotes the least common multiplicity. This result agrees with the above 
presented formulas, since x = gcd(K, A) = N/h with h = \cm{k,l). These double cosets 
are represented by ga = Clf with j = 0, 1, . . . , x — 1 if one subgroup is cyclic. If both 
subgroups are cyclic, then the next x representatives are gd = Uj, < j < x- The 
cardinality of double cosets just obtained does not depend on the representative g, so 
Eq. (P) can be simplified if one of subgroups is cyclic. To shorten the notation let us 
introduce two symbols: (i) the product Bl!/^{gd){U^gd\'H\Vi'E) is denoted as f{gd) {U, 
V, r etc. are omitted, but they can be easily determined from the context) and (ii) 
sd(A;, /) = lcm(A;, /)/ gcd(A;, /).[] Using these symbols one can write 

hc,,r,c,.sir) = sd{k,lfl'Y.fiCi,) + f{U'), (15) 

d=0 

V,.,c...(r) = /ic,,.,D,..(r) = (2sd(A:,/))V2 (16) 

The most complicated is the case when both subgroups U and V are dihedral, because 
double cosets obtained may have different lengths and there are three different values of 
the number of double cosets. In fact, we are not interested in their detailed structure, 
but we have to know the number of representatives g of left cosets gU contained in a 
double coset Vg^U. The appropriate rules can be derived if one notices that the action 
of V = Dk on the set of left cosets Djv/D^ can be performed in two steps: (i) the action 
of the cyclic subgroup Ck C considered above and (ii) the action of the generator Uq 
or Ui, depending on the type of V. The other elements Ug+p^, q = 0,1, < p < k, do 
not yield new results since t/g+p^ = UqCj^~^^'^ and the action of Cj^~^^'^ has been already 
taken into account. After the first step we have at hand x orbits represented by Cj^, 
< j < X- Each coset CfcC^D/ contains also X' = X/x representatives of left cosets gDi 
in the form C^^-' , < r < X'. From the multiplication rule (|^) we obtain (h = N/x)'- 

UoCif = ^U)^(h-i) or UoC^f = ^ ^Uy^(^h-i)+i ', 
UiCj^ = ''^'^U^{h-i) or UiC-'^ = ■'f^x('i-i)+i • 
The following four cases have to be discussed: 

• U = D°, V = D^: It follows from the first formula that cosets CfcC^D° and 
CkCpf~^'D^ are merged. There is one fixed point, namely j = 0. If x is even then 
there is another fixed point j = x/2. 

• U = Dj, V = Dl: In this case the last formula should be used so the results are the 
same as above. 

^The symbol 'sd' means the symmetric division and has been introduced by analogy to the symmetric 
difference of sets: sd(fc, I) is a product of primes, which appear in decomposition of either k or /; the 
authors have not found any special name for this symbol in available books on the number theory. 



Table 1: Representatives of cosets gU contained in the same double coset Vgll for the 
dihedral group Dg and U,V = D^, T)\, D2; underbraced elements belong to the same orbit 
under the action of the cychc subgroup C2 of V = D2 



u 

V 


D? D} D2 




{E}^ {E,CI} {E}^ 

{CI CD {CI CD 
{CD 


Dl 


{E,C,} {E} {E,Ce} 
{CI CD {Ce,CD {CD 

{CI CD {CI CD 
{CD 


D2 


{E,CD {E, CI CI CD {E} 



. U = Dj, V = DO; Orbits CfeC^D^ and CkCj^Dj are merged if f = -(1 + j) 
(mod x)- The fixed point is obtained for 2j = —1 (mod x), so it is possible for 
odd X only and in this case j = {x~ l)/2- 

• U = D^i, V = Y)\: In this case the congruence f = 1 — j (mod x) has to be 
satisfied to merge orbits and the fixed point is obtained for j = (% + l)/2, if x is 
odd. In this case it may be more convenient to use j = 1, 2, . . . , x/2 for x even and 
j = l,2,...,(x-l)/2,(x+l)/2forxodd. 

In Table |l] all possible cases are presented for = 6 and U,V = D^, Dj, and D2. Each 
entry is a set of representatives g of left cosets gU, which are contained in the same double 
coset VgdU . 

Taking into account the results presented above one obtains the following formulas for 
matrix elements: 

/ x/2-1 \ 

V>,DO..(r) = /^D>,Di..(r) = sd(fc,o^/M/(i?) + 2 f{c%) + f{ct)y,m 

x/2 

VM^,Di.s(r) = 2sd(A;,o'/'E/(c^); (is) 

d=l 

hY)0^r,T>kVs(^) = ^DiAtr,DO!/s(r) = ^D,Atr,Dfei/s (T) 

= sd(fc,/)i/2||/(E)+2^'|:^V(C^)j ; (19) 

x/2-1 

SVr,D0.s(r) = 2sd(A;,/)i/2 E f{C%)- (20) 

d=0 
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/ (x-3)/2 \ 

Sv,D...(r) = sd(A:,/)V2(2 Y: f{Cf,) + f{C^r'''^') \ ; (21) 
SMr,Di..(r) = sd(fc,/)^/^||2^'E^V(C'^^)+/(C'(,^+^)/^)j . (22) 

2.4 Transitive irreps 

Let us start with the cychc subgroups = (C^) with A = N/l. Since for all j = 
0, 1, 2, - 1 we have x^'''(C^) = 1 (see Eqs. (|,|10D), then 

n{Ai, R^''--^') = n{A2, = 1 . (23) 

In the case of even there is also another pair of one- dimensional irreps, namely Bi and 
i?2- It can be shown that both of them appears in the decomposition of R^^-Ci g^j-^^ only 
if A is even. The two-dimensional irreps Ep always appear twice, i.e. n{Ep, i?^^'*-'') = 2 if 
gcd(p, /) = I. Since p < N/2 then only subgroups with I < N/2 lead to decompositions 
containing Ep. 

The dihedral groups D/ have twice more elements then C; then the dimensions of 
R^N-'Di g^j,g times smaller. At first, one can notice that the irrep A2 does not appear 
in the appropriate decompositions. The case T = Ep is simply, since one obtains 

n{Ep, = n{Ep, i?°^^^0/2 = 1 if gcd(p, /) = / . 

The irreps Bi and B2 (admissible for even only) appears depending on the type of the 
group D^: 

n(Bi, i?°^^°?) = n{B2, i?°^^°') = 1 . 
The results obtained are illustrated in Table ^ for = 6. 

2.5 Reduction coefficients 

A general formula for the reduction coefficients, see Sec. 2.4 in yields that one- 
dimensional irreps F always lead to trivial reduction coefficients A^^ = 1. The indices 
7 = 1 and r = 1 can be omitted. So only the case T = Ep has to be considered in more 
details. The vector of irreducible basis related with the reduction Ep I U in this case can 
be written as 

\UEpr) = Al^-lEpx) + A^f^lEpy) , 

where r = 1 for the dihedral (D^) and r = 1, 2 for the cyclic (C/) subgroups U, respectively. 
In the first case we have to distinguish subgroups D° and D^^. The standard choice of the 
basis {\Epx), \Epy)} and Eq. ([T^) yield the matrices: 

/I 0\ p(Tj^( cos{2np/N) sin(27rp/Ar) \ 

I -1 r P^^' [ sin{27ip/N) -cos{2np/N) ) ' 



Table 2: Decompositions of transitive representations into irreps for Dg (subgroups are 
determined by their generators) 



u 


















3 El ^9 © 2i?i © 2E'} 


D? = 


(Un) 




^ Bli 


B £^1 © 


Di = 




Ai 6 




3 ^1 © ^2 


C2 = 


Ai 6 




3 2^2 


D2 = 


f/o) 


v4i 6 


3^2 




C3 = 




3^2 6 


3 fil © ^2 


D^ = 


Ai 6 


3fii 




DH 


Ai 6 


3^2 




C6 = 


(Ce) 


Ai 6 


3^2 




D6 = 


{CM) 









Therefore, one obtains (for representations appearing in the decompositions of 
i.e. for gcd(p, V) = I): 

Ax''''' = l, Ay°''' = 0; (24) 

= cos((/./2) , A?' '''' = sm{^/2) , (25) 

where = 2ttp/N = 2ttpi/\ with pi = p/l. The cychc subgroups yield r = 1, 2 and it can 
be shown that there are no restrictions imposed on the coefficients A^}.^p, as long as they 
determine orthonormal vectors \CiEpl), \CiEp2). The simplest choice is to put 

A':1''- = 1, A^l'''' = 0; (26) 

A^^^'' = 0, Ag^- = 1. (27) 

Now we have to calculate the coefficients B^^^{gj) for the representatives of left cosets 
QjU. Let us start from one-dimensional irreps F = Ai, A2, Bi, B2. Since the coefficients 
A^^ = 1 and r,s = 1, then 

B'^'^'ig,) = T{g,) = x^{g,) . (28) 

Characters of these representations are given by Eqs. (^-(0). In particular, the unit 
irrep Ai always yields B^^^^{gj) = 1. 

The two-dimensional irreps Ep appear in decompositions of R^n-Ci g^^^j ji^^n-'Di ^^j. 
I < N/2 and p = Ipf, if U or V is the cyclic subgroup then Ep appears twice. Moreover, 
there are two types of the dihedral subgroups for even. It leads to 16 different cases of 
coefficients B^^ ^p, which will be presented below. We are interested in the coefficients 
B^^ when the irrep Ep appears in the decompositions of R^^-U g^j^^ R^^-V _ ^^^^ 

1 o 



in Sec. 2^, V = Ck or and U = C/ or D^. Then the index p of the irrep Ep has to 
satisfy both conditions: p = pkk and p = pil. Therefore, p = lcm{k,l)p' = hp' = Np'/x 
and if = 2ttp/N = lup' jx. 

Let subgroups f/, V be cychc. Then both repetition indices r, s = 1, 2, so there are four 
different coefficients B^^ ^^icj^). We have chosen elements and Uj, j = 0, 1, . . . , A — 1 
as representatives of the left cosets g^Qi with representation matrices 

EAC'^) = ( ] , EM) = ( '"^"'-^ ] ■ 

^ \ sm jip cos j(p I F\ jj \ sm jip — cos j(p I 

Taking all these relations into account one obtains 

B^f'''^^{Cj,)=B^''''''''{Cj,)=B^f'''^{U,) = -B^f''''''{U,) = cosj^, (29) 
B^f'^^^iCj,) = -B^f'^'^^iC],) = Bff^'^-iU,) = B^f'^'^iU,) = sinjyp . (30) 

When the subgroup V is dihedral than the repetition index can be omitted, but two types 
of dihedral subgroups have to be considered. It leads to the following formulas: 

5ii?°''^(C^)=5^i?°''^(f/,) = cosj^, (31) 

5^:?''^(f/,) = -5S°''^(C'^) = sinjy,, (32) 

B^lT'^'^Cj,) = B^f,'^''%U,) = cos(j-l/2)^, (33) 

BrLT'iU,) = -B'^LT'iCj,) = sin(j-l/2)y,. (34) 

The case U = D;, V = is a bit easier to discuss, since the left cosets gjDi are represented 
by the elements Cj^ only, 

B^l^'^'^iCj,) = cosjV, (35) 
B^l^^'^igj) = sin jip, (36) 

5fJ^^''(C^) = cos(j + l/2)¥, , (37) 
i?fi'=^''(^7,) = sin(j + l/2)^. (38) 

In the last four cases both subgroups are dihedral, so gj = Cj^ and different types of 
subgroups are considered. The final results follows: 

5D?D2^.(C^) = cosjV, (39) 

i?D^'°°^-(C^) = cos(j + l/2)(^ , (40) 

B^'^l^^iC'^) = cos(j - l/2)y, , (41) 

B^^^'^^^iCj,) = cos Jip. (42) 

The equation presented together with Eqs. (|T5|-^) enable us to determine matrix 
elements for any operator commuting with P{g). 



Table 3: Dimensions of subspaces of given magnetization M > and non-ordered parti- 
tions related with them 



M non-ordered partitions [/c] \^[k]\ dimL^ 

~0 [303], [222], [141], [060] 20, 90, 30, 1 lIT 

1 [213], [132], [051] 60, 60, 6 126 

2 [204], [123], [042] 15, 60, 15 90 

3 [114], [033] 30, 20 50 

4 [105], [024] 6, 15 21 

5 [015] 6 6 

6 [006] 1 1 



3 Example 

Let us consider a ring of six spins s = 1 with the symmetry group Dg and the antiferro- 
magnetic interactions of the nearest-neighbors (J = —1) 

6 

'H = ^Sj- Sj+i , 6 + j = j ; 
i=i 

The dimension of the space of states is (2s + 1)^ = 729, so it is not so large, but it is 
sufficient to illustrate application of the method presented in this paper. 

Since in this case the magnetization is a good quantum number, then at first we 
consider ordered and non-ordered partitions of = 6 into no more than 2s + 1 = 3 
parts. There are seven ordered partitions: [600], [510], [420], [411], [330], [321], and 
[222]. They lead to 28 non-ordered partitions, which were presented in the previous paper 
1^ in Table 1. On the other hand, some simple combinatorial considerations yield the 
dimensions of subspaces with given magnetization M |P, ITof. The results obtained are 



collected in Table ^ Due to the time-reversal symmetry it is enough to consider the states 
with M > 0. 

The non-ordered partitions determine orbits of the symmetric group Eg. The action 



of Dg C Sg decomposes them into orbits labeled by subgroups [/ C Dg (see Sec. ^ 
and Fig. ||). These decompositions depend only on the type of [fc], i.e. on the ordered 
partition [k], and are presented below (the number before a subgroup symbol U denotes 
the multiplicity of orbits of the same type U): 

[6]:D6; [51]: D?; [42] : D2, D?, Dj; [33] : D^, D?, Ci; 
[411] :D°,2Ci; [321] : 2D°, 4Ci; [222] : C2, 3D?, 3Di, 4Ci . 

These results yield classification of orbits by the magnetization M and the type U : 

M = 6:D6; M = 5 : D? ; M = 4 : D2, 2D?, D} ; M = 3 : D°, 2D?, 3Ci ; 



Table 4: Classification of states by the magnetization M > and the irrep T 



M 


Ai 


^2 


Bi 


B2 


El 


E2 


dim. Lm 


u 


1 


n 


n 


n 


n 


n 


1 

L 


5 


1 





1 





1 


1 


6 


4 


4 





2 


1 


4 


3 


21 


3 


6 


3 


6 


3 


8 


8 


50 


2 


12 


4 


8 


6 


16 


14 


90 


1 


13 


8 


13 


8 


21 


21 


126 





18 


8 


13 


10 


22 


24 


141 


Total 


92 


38 


73 


46 


122 


118 





M = 2: 2D2, 4D°, 2D}, 4Ci ; M = 1 : 5D°, 8C1 ; M = : Dg, D°, C2, 5D?, 3D}, 7Ci . 

Now, taking into account the decompositions presented in Table ^, we can write down 
the multiplicities n(r, P I Lm) for each M > 0. The numbers presented in Table ^ give 
us the dimension of the eigenproblem 'H\ip) = E\iIj) restricted to the states with a given 
magnetization M and the symmetry described by F. The largest dimension is reached 
for M = and T = E2. However, when we are interested in the ground state only, 
then r = Ai and the largest eigenproblem has dimension 18 ~ 729/40. We can decrease 
these numbers even more taking into the square of the total spin S^. Since eigenvalues 
of this operator are S{S + 1) for < S* < 6, then it is enough to find the corresponding 
eigenvectors and express the Hamiltonian matrix in this new basis. Table ^ provides us 
also with dimensions dimLM=s,r of subspaces with a given total spin S and the maximum 
magnetization M = S corresponding to an irrep F. It is enough to subtract adjoining 
rows of this table (except for the first one since dim Lm=s=Ns = dim Lm=ns), because 
sates with M = S are these with a given M except for those with S = M + 1, which are 
listed in the previous row. For example, the dimensions dimL^^^^^^^ are 5, 1, 6, 2, 3, 0, 
and 1 for 5 = 0, 1, . . . , 6, respectively. In the same way one obtains that the subspace of 
states with M = S = has the dimension dimLM=5=o = 15 and 

P i Lm=s=o = 5v4i © 2B2 © ^1 © 3^2 , 

so the dimension of the Hamiltonian eigenproblem in this case is at most five. The overall 
maximum is reached for M = S = 2 and T = Ei and equals eight, i.e. about one ninetieth 
of the original dimension dimL = 729. 

Let us restrict the further considerations to the case M = S = and T = Ai. There 
are 18 orbits grouped into four classes labeled by the non-ordered partitions [060], [141], 
[222], and [303]; their representatives are presented in Fig. |. Using the symbols \U[k]^), 
where U is a stabilizer, [k] — a non-ordered partition, and fi — a configuration expressed 
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Figure 3: Orbit representatives for iV = 6, s = 1 and M = 0; full and empty circles 
denote = +1, — 1, respectively 

as a series (mi, . . . with ± standing for ±1, respectively, these representatives can 
be written as: 



|1) 
1 3) 
|5) 
1 7) 
|9) 

111) 
1 13) 

|15) 

1 17) 



IDs [060 
|Ci [141 
IC2 [222 
|D; [222 
|D} [222 
|D} [222 
jCi [222 
|Ci [222 
|D° [303 



(0,0,0,0,0,0)), 
(+,-,0,0,0,0)), 
(+,0,-,+,0,-)) 
(0,+,-,0,-,+)) 
(+,+,0,-,-,0)) 
(-,-,+,0,0,+)) 
(+,0,+,0,-,-)) 
(+,+,0,0,-,-)) 

(+,+,-,-,-,+)), 



|2) 
|4) 
|6) 
|8) 
10) 
12) 
14) 
16) 
18) 



ID? [141] 
|Ci [141] 
ID? [222] 
|D? [222] 
|D} [222] 
jCi [222] 
jCi [222] 
ID3 [303] 
|Ci [303] 



+,0,0, -,0,0)), 
+,0,-, 0,0,0)), 
+,0,-,+,-,0)) 
-,+,0,-,0,+)) 
+,+,-,0,0,-)) 
+,+,0,-,0,-)) 
+,-,+,0,0,-)) 

+,-,+,-,+,-)), 
+,+,-,+,-,-)). 



Each orbit gives rise to only one state |j; A\) = |/), j — 1,2, ... , 18, with symmetry F = 
Ai constructed as a normahzed sum of all |D6|/|C/| states in a given orbit 0{U, [A;], /x) = Oj. 

Since the Hamiltonian considered is isotropic, then we can use the total spin number 
S to label eigenstates. The operator can be expressed as 



6 5 6 6 3 

S' = Esi+2E E Si-Si, ^12 + 2n + 2j2si-Si+2 + 2Y,Si 

i=l i=l i'=i-\-l i=l i=l 

The diagonal parts of 7i and are given as 



Sj+3 ■ 



H 



diag 



g2 

diag 



E«^-«' 



22 



i=l 



12 + 27idiag + 2 E + 2 E 

i=l i=l 
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In the first case (for TYdiag) one obtains the following 18 eigenvalues: 0, 0, -1, 0, -2, -2, 
-2, -2, 2, -1, -1, 0, 0, -3, 1, -6, 2, and -2. The case of operator is even more simple, 
since the diagonal values depend only on the partition [k] and one obtains four numbers: 
12, 10 (orbits 2, 3, 4), 8 (orbits 5 through 15), and 6 (orbits 16, 17, 18). 

To consider off-diagonal terms containing step operators sf we start with a graph of 
partitions: 

D_0 

[060] [141] [222] [303] 

This graph leads a block structure of off-diagonal parts of the Hamiltonian and (sym- 
metric) matrices: 

/ C \ 
A V 
B £ ' 

V £^ J 

where rows and columns are labeled by non-ordered partitions and f denotes the Hermi- 
tian conjugation. All these blocks may contain non-zero elements for F — Ai, because this 
representation appears in the decomposition of any transitive irreps. More detailed struc- 
ture can be obtained if orbit stabilizers are taken into account (now rows and columns 
are labeled by pairs {U, [k]); borders of the blocks introduced above are denoted by hnes): 








C2 








cl 


A, 


As 


Vl 


V2 


Vs 


V, 











4 


A2 


V5 


Ve 


Vr 


Vs 








vl 


Vl 


Bi 


B, 


Be 


B7 


£1 


£2 


S3 





vl 


vl 


Bl 


B2 


Bs 


B, 


£a 


£, 


s. 


vl 


v\ 


Bl 


Bl 


B3 


Bio 


£1 


£s 


£^ 






^8 


Bl 


Bl 


Bio 




£10 


£11 


£12 








4 


Si 




£^ 
^10 
















Si 

si 


Si 
Si 


F^ 

F^ 


F^ 
^11 

F^ 
^12 










Depending on the representation considered some blocks have to contain only zeroes. For 
example, when V = A2, then only blocks A2, Bi, B4, Bj, P5, Pg, £3, £12 may contain 
non-zero entries. It follows from the fact that the irrep A2 appears in the decomposition 
of a transitive representation R^'^'-'^ if f/ is a cyclic subgroup of Djy. 

The classification of orbits with M = by type U yields that there are six different 
types of orbits, namely Dg, D3, C2, Dj, D{, Ci, so there are 21 different double cosets to 
be considered.. However, the matrix structure presented above excludes cosets DqqDq, 
DegDl DegD\, Dg^Cs, and D^gDl 

Two notes are in place: (i) since we consider the case F = Ai, so all coefficients 
^uvr _ "1^. ^--^ ^Yie off-diagonal terms contains products sfs^, which acting on kets 
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|mjmj/) gives zero or 2|mj + 1 rrij/ — 1), therefore terms {sfs^ + s~s^)/2 are replaced by 
cr^o-^ + '^i'^t^ where af^\ ± 1) = |) and erf |mj) = |mj ± 1), otherwise. 

Let us consider, as an example, the 16th column of an operator matrix. This column 
corresponds to the two-element orbit containing the Neel configurations |+,— ,+,—,+,— ) 
and | — , +, — , +, — , +) with a stabilizer D3 and labeled by the non-ordered partition [303]. 
Therefore, non-zero entries in this column, for operators constructed of bilinear terms, 
can appear only in the block S, i.e. rows labeled by the partition [222] (except for the 
diagonal terms Z]j,j' -sf sf/). So, one has to consider four sets of double cosets DI\Dq/U 
with U = C2, D5, and Ci. In all cases k = 3 and k = 3; in the first case {U = C2) 
/ = 2, A = 3, so sd(A;,/) = 6 and x = whereas in the other cases we have / = 1, 
A = 6, sd{k, /) = 3 and x = 2. Eqs. (p!6|), (p!7D, and ( |20D yield the following formulas (the 
repetition indices r, s may be omitted): 

hc,,,DoMi) = 2^f{E); hc,,,noMi) = V6{f{E) + f{C,)); 
hBl,,DoM^) = 2v^/(E); V>,DO.(Ai) = V3{fiE) + f{Ce)). 

Therefore, one has to take into account the states for 5 < j < 15 and P{CQ)\j) for 
j = 6, 7, 8, 12, 13, 14, and 15. Only two of them are obtained from |16) under the action 
of bilinear terms crf (7^7, namely 

PiCe)\7) = |+, 0, +, -, 0, -) = a+a^\+, -, +, -, +, -) , 

and 

|14) = |+,-,+,0,0,-) =a+cr5-|+, -,+,-,+,-). 
Hence, the non-zero off-diagonal terms in the sixteenth columns are 

7^14,16 = ; si,, = 2V3 , s?4^i6 = 2 . 

All other elements of operator matrices can be obtained using the following formulas: 

/iCiM,D6.(^i) = Vm.D6.(^) = 2v^/(E); 

2 

d=0 
d=0 

/^DV,DO.(^i) = /^DiM.Di.(^i) = f{E) + f{Cl) + 2{f{Ce) + f{Cl)); 

d=0 

d=Q 

1 o 



Table 5: The Hamiltonian (the upper triangle) and (the lower triangle) matrices for 
N = 6 spins s = 1 in the basis {\j') = \j',Ai) | 1 < j < 18} (the lines are borders of 
blocks A, B, C, V, and £) 









2V2 




















\/2 \/2 








-1 


2 






V2 






^2 V2 






2 








6 










V2 












1 


1 










V6 


5 






-2 


2 




2 






















2v^ 


7 






-2 
















V2 








2^3 


2^2 


4 


7 






-2 










72 V2 














V2 


V2 


4 






-2 




























2 


4 






2 


















2 










4 






-1 














x/2 






72 


V2 


2 






4 






-1 












x/2 




2 




V2 










4 












1 












V2 


2 






2 




4 










1 










V2 


V2 


2 


2 










4 






-3 








1 




V2 


1 


2 










V2 




5 






1 


■ 


1 




x/2 


1 


2 
















5 






-6 








V2 


3 






V2 2V2 V2 








1 


1 


6 




2 








1 


2 










2^2 72 V2 


2 


2 




4 




-2 
































3 


















1 




2 






V2 






■ 3 




















V2 


1 


1 


1 


1 




3 



The matrices of the Hamiltonian Ti and the square of total spin are presented in 
Table 

The irrep Ai has not been chosen merely due to the simplest form of the coefficients 
^c/yAi^^^^ but in order to illustrate the other way of determination of matrix elements 
in this case. Let an orbit Oj consists of Uj = \G\/\Uj\ elements [Uj is a stabilizer of the 
configuration \j) = \Uj [k] /i)). Hence, the fully symmetric state = |j; Tq) is given as 

\j') = -=T.Pi9i)\j), (43) 

where gt are representatives of left cosets QiUj G G/Uj. Since a stabilizer of each element 
in the orbit Oj has \Uj\ elements than 



I/) = WV^ ^ ^^3^)\J) = ^1: Pi9)\j) ■ 

\'-'j\V i g&G \^\ g€G 

Let Ti.' denote the off-diagonal part of a Hamiltonian Ti, i.e. Ti' = Ti — Tidiag- The 



symmetry operators P{g) commute with a Hamiltonian Ti.', so 




E Pi9)n'\j) . 



G 



geG 



To simplify the further considerations we hmit ourselves to the cases s = 1/2,1, when 
sfs^\mi . . . , mj/) = \mi + 1, m^/ — 1) or cr^a^\mi, rriii ■ ■ ■) = Ifrii + 1, mj/ — 1), since in the 
case of larger spin numbers more detailed discussion is necessary. Under this assumption, 
each of off-diagonal terms of the bilinear Hamiltonian gives the empty sate |) or one 
configuration without any coefficients. Let denote a number of configurations obtained 
from I j) (by the action of H') and belonging to the orbit Ok- Then J2g P{g)'H'\j) contains 
\G\njk configurations from the orbit Ok- Since the state \k') = \k; Ai) are constructed 
from all Uk elements of the orbit O^, then \G\njk has to be divisible by and each element 
the orbit Ok appears \G\njk/nk times in the considered sum X)g P{9)'H'\j)- Therefore, one 
can construct \G\njk/ ^/nk vectors \k'). As the result we obtain 



In the example presented above \j) = |+, —,+,—,+, —) and the off-diagonal terms of the 
Hamiltonian give us six configurations belonging to the orbit O14 (with ni4 = 12). The 
additional terms in operator lead to three configurations in the 6-element orbit O7. 
Therefore, 



what agrees with the previous results. 

The further considerations go beyond the scope of this paper, since there are some 
numerical rather than combinatorial problems. However, we are very close to the solution 
of the eigenproblem and it is interesting to calculate the ground state energy. In the 
case of isotropic Hamiltonians, when S* is a good quantum number, we determine an 
orthonormal basis for each space Lm=s,Ai solving the eigenproblem S^|0) = S{S + 1)|0) 
for 5 = 0, 1, 2, 3, 4, 6 (S* = 5 may be omitted since there are no states with M = S = 5 
and T = Ai). It can be done, for example, using such packages as MATHEMATICA® 
or numerical libraries for calculations with an arbitrary precision as GMP Then, 
we transform the Hamiltonian to this new basis, what leads to a quasi-diagonal matrix 
consisting of blocks with dimensions 5, 1, 6, 2, 3, 1 for S" = 0,1,2,3,4,6, respectively. 
It is easy to solve the eigenproblems for such small blocks and, for example, in the case 
S = we obtain five eigenvalues: -8.6174, -4.7979, -3.3302, -0.3391, 3.0846, where the 
first one corresponds to the ground state with energy per spin -1.43624. In the case of 





such small eigenproblems one can also determine eigenstate^ and, therefore, it is possible 
to calculate spin correlations in the ground state. 

4 Final remarks 

In this paper we have studied properties of the dihedral groups D^f to obtain analytical 
formulas for number of double cosets and their elements. These results are necessary 
to exploit in full symmetry properties of a given spin operator and to obtain its matrix 
elements in the symmetry adapted basis. The groups considered are important in inves- 
tigation of magnetic macromolecules forming rings or other structures with isomorphic 
symmetry group {e.g. crown-like structures). In a sense this paper, together with the 
preceding one [^, continues works done for the cyclic groups by T. Lulek and one of the 
authors (WF) . In those articles the arithmetic structure of a number of nodes was 
mostly taken into account and some combinatorial object were considered. The present 
series of papers is mainly devoted to combinatorial and group-theoretical objects and the 
most important result is the formula for matrix elements expressed as the sum over repre- 
sentatives of double cosets. In the case of dihedral groups the approach proposed has two 
advantages: (i) all formulas can be obtained in analytical form and they can be used in 
computer programmes; (ii) in one dimension a loop over double coset representatives is, 
in general, shorter than a loop over all off-diagonal elements of bilinear Heisenberg Hamil- 
tonian, what has been discussed in the previous paper 0; moreover, it is easy to modify 
Hamiltonian introducing, for example, interactions with the next nearest neighbors or 
alternating exchange integrals. 

It should be underlined, however, that the final diagonalization of a given Hamiltonian 
is (almost) always done numerically, with all round-off errors caused by the float-point 
representation of numbers. Secondly, the method proposed only slightly increase the size 
of system which can be investigated. The so-called 'combinatorial explosion' leads to very 
large eigenproblems for quite small systems: e.g. a system consisting of twelve spins s = 2 
yields more than 240 millions Ising configurations, so there are about ten millions orbits 
of the dihedral group D12; more than 800 thousands orbits correspond to M = 0, so this is 
the approximate of eigenproblem dimension. On the other hand, not so much smaller 
system with = 8 leads to about 2,500 orbits with M = 0, what can be investigated using 
modern computers. A number in the same range is obtained for A^ = 19 and s = 1/2, 

^Due to its form the eigenproblem for can be solved exactly in many cases. This leads to an exact 
form of blocks in the Hamiltonian matrix. However, the Hamiltonian eigenproblem, except for some 
special cases, has to be solved numerically with the precision depending on software and hardware used. 
Therefore, all numbers obtained give only approximate values of energies, correlations, vector coefhcients 
etc. 
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but we hope to determine the ground state energy for the Heisenberg antiferromagnetic 
ring with N = 24 nodes, where a system of about 50 thousands hnear equations has to be 
solved. However, if one would like to determine the thermodynamic properties then all 
states have to be taken account, among others these with the symmetry F = Ep, which 
are double-fold degenerated, so the dimension of eigenproblem is approximately two 
times larger than that for one- dimensional irreps. 
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Figure captions 

Figure 1: Numbering of ring nodes and elements of the group Dg (the six-fold axis Cq 
is denoted by the full hexagon) 

Figure 2: Lattice of conjugated subgroups for the dihedral group Dg; lines correspond 
to the group-subgroup relations 

Figure 3: Orbit representatives for = 6, s = 1 and M — 0; full and empty circles 
denote rrij — +1, —1, respectively 
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